Self-induced tunable transparency in layered superconductors 
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We predict a novel nonlinear electromagnetic phenomenon in layered superconducting slabs irra- 
diated from one side by an electromagnetic plane wave. We show that the reflectance and transmit- 
tance of the slab can vary over a wide range, from nearly zero to one, when changing the incident 
wave amplitude. Thus changing the amplitude of the incident wave can induce either the total 
transmission or reflection of the incident wave. In addition, the dependence of the superconductor 
transmittance on the incident wave amplitude has an unusual hysteretic behavior with jumps. This 
remarkable nonlinear effect (self-induced transparency) can be observed even at small amplitudes, 
when the wave frequency lu is close to the Josephson plasma frequency ojj. 

PACS numbers: 74.78. Fk, 74.50.-|-r, 74.72.-h 
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I. INTRODUCTION. 

The recent growing interest in the unusual electrody- 
naniic properties of layered superconductors (see, e.g., 
the very recent review Ref. [l|) is due to their possi- 
ble applications, including terahertz imaging, astron- 
omy, spectroscopy, chemical and biological identification. 
The experiments for the c-axis conductivity in high- Tc 
Bi2Sr2CaCu2 08-1-5 single crystals justify the use of a 
model in which the very thin superconducting Cu02 lay- 
ers are coupled by the intrinsic Josephson effect through 
the thicker dielectric layersr^—. Thus, a very specific type 
of plasma (the so-called, Josephson plasma) is formed in 
layered superconductors. The current-carrying capabil- 
ity of this plasma is strongly anisotropic, not only in the 
absolute values of the current density. Even the physical 
nature of the currents along and across layers is quite dif- 
ferent: the current along the layers is the same as in usual 
bulk superconductors, whereas the current across the lay- 
ers originates from the Josephson effect. This Josephson 
current flowing along the c-axis couples with the elec- 
tromagnetic field inside the insulating dielectric layers, 
causing a specific kind of elementary excitations called 
Josephson plasma waves (JPWs) (see, e.g., Ref. [H and 
references therein). Therefore, the layered structure of 
superconductors favors the propagation of electromag- 
netic waves through the layers. 

The electrodynamics of layered superconductors is 
described by nonlinear coupled sine-Gordon equa- 
tionsii^"— . This nonlinearity originates from the nonlin- 
ear relation J oc simp between the Josephson interlayer 
current J and the gauge-invariant interlayer phase differ- 
ence if of the order parameter. As a result of the nonlin- 
earity, a number of nontrivial nonlinear phenomenapii"— 
(such as slowing down of light, self-focusing effects, the 
pumping of weaker waves by stronger ones, etc.) have 
been predicted for layered superconductors. The nonlin- 
earity plays a crucial role in the JPWs propagation, even 
for small wave amplitudes, \(p\ ^ 1, when simp can be 



expanded into a series as (p — p^ /6), for frequencies close 
to the Josephson plasma frequency. 

In this paper, we predict a novel and unusual strongly 
nonlinear phenomenon. The reflectance and transmit- 
tance of a superconducting slab being exposed to an in- 
cident wave from one of its sides depend not only on the 
wave frequency and the incident angle, but also on the 
wave amplitude. If the frequency a; of irradiation is close 
to the Josephson plasma frequency wj, the transmittance 
of the slab can vary over a wide range, from nearly zero to 
one. Therefore, a slab of fixed thickness can be absolutely 
transparent (when neglecting dissipation) for waves of 
definite amplitudes, and nearly totally reflecting for other 
amplitudes. This unusual property can be described as 
self-induced tunable transparency. Moreover, the depen- 
dence of the transmittance on the wave amplitude shows 
hysteretic behavior with jumps. Tunable clcctromagnct- 
ically induced transparency is also now being studied in 
various contexts, including superconducting circuits and 
other lambda-type atoms (see, e.g., Refs. Hal 171 1. 

The large sensitivity of the slab transmittance to the 
wave amplitude can be explained using a very clear phys- 
ical consideration. Let us consider a wave frequency uj 
which is slightly smaller than the Josephson plasma fre- 
quency ujj. In this case, linear JPWs cannotpropagate 
in the layered superconductor (see, e.g., Ref. [l|). So, the 
slab is opaque for waves with small amplitudes, and the 
transmittance is exponentially small due to the skin ef- 
fect. According to Refs. Illl - USl the nonlinearity results in 
an effective decrease of the Josephson plasma frequency 
and, thus, nonlinear JPWs with high-enough amplitudes 
can propagate in the superconductor. Moreover, chang- 
ing the amplitude of the incident wave one can attain 
the conditions when the slab thickness equals an integer 
number of half- wavelengths. Under such conditions, the 
slab becomes transparent, with its transmittance equal 
to one. 

The paper is organized as follows. In section II, we dis- 
cuss the geometry of the problem and present the main 
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B. Electromagnetic field in the vacuum 

The magnetic field H" in the upper vacuum semispace 
(z > D) can be represented as a sum of the incident and 
reflected waves with amplitudes Hg and Hr , respectively. 
The field H^ in the vacuum semispace below the sample 
{z < 0) is the transmitted wave with amplitude Ht . The 
upper (ff") and lower (ij') fields and can be written in 
the following form, 

iJ" = Ho cos[fea:a; — Lot — kz{z ~ D)] 



FIG. 1: (Color online) Geometry of the problem. A slab of 
layered superconductor is irradiated from the upper side by a 
p-polarized electromagnetic wave. 



Hr cos[A:2;a:; — ujt + kz{z — D) + a], (2) 



7j' = Ht cos{kxX — (jjt — kzZ + /?). 



(3) 



equations for the electromagnetic fields both in the vac- 
uum and in the slab of layered superconductor. In sec- 
tion III, we express the transmittance T in terms of the 
amplitude of the incident wave and analyze this depen- 
dence in two cases: when the frequency u of the inci- 
dent wave is either larger or smaller than the Josephson 
plasma frequency loj. In both cases, we study the un- 
usual hysteretic features of this dependence. The results 
of numerical simulations support our theoretical predic- 
tions. 



Here 
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are the components of the wave- vector k; of the incident 
wave, a and /3 are the phase shifts of the reflected and 
transmitted waves. Using Maxwell equations one can 
derive the electric fleld components in the vacuum: 

£■" = —Hocos6cos[kxx — ojt — kz{z — D)] 



II. SPATIAL DISTRIBUTION OF THE 
ELECTROMAGNETIC FIELD 

A. Geometry of the problem 

We study a slab of layered superconductor of thick- 
ness D (see Fig. [T|) . Superconducting layers of thickness 
s are interlayed by insulators of much larger thickness d 
{s <^ d). We assume the number of layers to be large, 
allowing the use of the continuum limit, and not to con- 
sider the spatial distribution of the electromagnetic fleld 
inside each layer. The coordinate system is chosen in 
such a way that the crystallographic ab-plane coincides 
with the xy-plane, and the c-axis is along the z-axis. The 
plane z = corresponds to the lower surface of the slab. 

A monochromatic electromagnetic plane wave of fre- 
quency Lu is incident on the upper surface of the slab, 
and it is partly reflected and partly transmitted through 
the slab. We consider the incident wave of the transverse 
magnetic polarization, when the magnetic field is parallel 
to the surface of the slab, 



Hr cos 9 cos[kxX — Ljt + kz{z — D) + a], (5) 



E^ ~ —Ho sin 6 cos[kxX — ut — kz{z — D)] 



Hr sin 9 cos[kxX — ut + kz{z — D) + a], (6) 



S^ = -HtCos9cos{kxX-ujt-kzZ + f3), (7) 



eI = - Ht sin 9 cosik.j:X - Lot - kzZ + fi). (8) 



C. Electromagnetic field in the layered 
superconductor 

The electromagnetic field inside a layered superconduc- 
tor slab is determined by the distribution of the gauge- 
invariant phase difference ip{x,z,t) of the order parame- 
ter between the layers (see, e.g., Ref. [l|). 



E = {E,,0,E,}, H = {0,H,0}. 



(1) 
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The incident angle 9 is considered to be not close to zero, 
so that both components k^ and kz of the wave- vector k^ 
are of the order of lo/c. 
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Kbd'H^ 
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Here Ho = $o/27rdAc, $o = ncH/e is the magnetic flux 
quantum, Xab and Ac = c/ujje^^^ are the London pene- 
tration depths across and along the layers, respectively. 
The Josephson plasma frequency is defined as 



wj 



SnedJc 



he 



(10) 



where Jc is the critical value of the Josephson current 
density, e is the permittivity of the dielectric layers in 
the slab. We omit the relaxation terms because, at low 
temperatures, they do not play an essential role in the 
phenomena considered here. 

The phase difference (p obeys a set of coupled sine- 
Gordon equations, that, in the continuous limit, takes 
on the following form (see, e.g., Ref. [l| and references 
therein) : 
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0. (11) 



In this paper, we study the case of weak nonlinear- 
ity: when the Josephson current density Jc sin ip can be 
expanded as a series for small ip, up to the third order, 
Jc simp « Jc{(p — V'^/6). We consider frequencies ut close 
to ujj and introduce a dimensionless frequency. 
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close to one. In this case, in spite of the weakness of the 
nonlinearity in Eq. (|lip . the linear terms nearly cancel 
each other, and the term Lp'^ plays a crucial role in this 
problem. Moreover, when the frequency ui is close to 
the Josephson plasma frequency wj, one can neglect the 
generation of higher harmonicai^^. 

It should be also noted that the nonlinearity provides 
an effective decrease of wj. Indeed, the expression in the 
square brackets in Eq. (jlip can be presented in the form 
[{ujj^)~^d'^(p/dt'^ + ip], where 
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For not very small ip, the frequency of the incident 
wave can be greater than the effective Josephson plasma 
frequency Wj and, therefore, the nonlinear Josephson 
plasma waves can propagate across the superconducting 
layers. 

The z-component of the electric field induces a charge 
in the superconducting layers when the charge compress- 
ibility is finite. This results in an additional interlayer 
coupling (so-called, capacitive coupling). Such a cou- 
pling significantly affects the properties of the longitu- 
dinal Josephson plasma waves with wave-vectors per- 
pendicular to the layers. The dispersion equation for 
linear Josephson plasma waves with arbitrary direction 
of the wave-vectors, taking into account capacitive cou- 
pling, was derived in Ref. [l^ According to this dis- 
persion equation, the capacitive coupling can be safely 



neglected in our case, when the wave-vector has a com- 
ponent kx '^ co/c along the layers, due to the smallness 
of the parameter a = Rj^e/sd. Here Rd is the Debye 
length for a charge in a superconductor. 

We seek a solution of Eq. pT|) in the form of a wave 
running along the x-axis, 

ip{x,z,t) = a{z)\l - n'^l^^^ sm[kxX - ijt + r]{z)], (12) 

with the z-dependent amplitude a and phase 77. We in- 
troduce the dimensionless z-coordinate, 



KZ ^cf^x 



(13) 



and the normalized thickness of the sample S = Dk/ Xab- 
Substituting the phase difference ip given by Eq. (J12p 
into Eq. (|lip . one obtains two differential equations for 
the functions 77(C) and a{(). The first of them is 



^'(0 



L 



hHO 



(14) 



where L is an integration constant, prime denotes deriva- 
tion over C, and 



HO = ±a(C) 



a'iO 



(15) 



The sign in this equation is plus for il < 1 and minus 
for ri > 1, i.e., it is opposite to the sign of the following 
important parameter, the 

frequency detuning = fl — 1. 

The coupled sine- Gordon equations PTj) give also the 
differential relation for /i(C): 



1^ 



h 



(16) 



Equations ((H]), dH]), dH]), and ^^ allow one to calcu- 
late the distribution of the phase difference ip{x, z, t) and 
then, using Eq. ©, the electromagnetic field inside the 
superconducting slab. 



III. TRANSMITTANCE OF THE 
SUPERCONDUCTING SLAB 

A. Main equations 

In this section, we analyze the transmittance T of a 
slab of layered superconductor. We rewrite the expres- 
sions for the magnetic field H^ and for the cc-component 
E^ of the electric field inside the slab using Eqs. © and 
Eq. (USD, 

H'ix, C, t) - ^Ho^l—^h{C) cosikxx -Lut + 77(C)) , 



|l-172|cos6' 



E:^ix,c,t) = nor^ 



X [h{C) sin [kxX - ujt + r]{C))]' . 
Here we introduce the parameter 



(17) 



Ac\/ecos0' 



which is usually small for layered superconductors. Now 
we can find the unknown amplitudes of the reflected and 
transmitted waves by matching the magnetic fields and 
the cc-components of the electric field at both interfaces 
(at z — and z ~ D) between the vacuum and the 
layered-superconductor. Using Eqs. p7|) for the fields in 
the superconductor and Eqs. 0, (j31), ©, and ([7]) for the 
fields in the vacuum, we obtain the following three equa- 
tions for the amplitudes a(0), a{S) and their derivatives 
on both surfaces of the layered superconductor: 



HS) 



h[5)) 

a'(0) = 0. 



+ r'[h'iS)]' = 4hi 



Here 



ho = 



Hn 



Ha |l-f^2 



(18) 

(19) 
(20) 



(21) 



is the normalized amplitude of the incident wave. These 
three equations, together with the coupled sine- Gordon 
equations (|16p. determine the integration constant L for 
each amplitude of the incident wave h^ (see Eq. (fT4|) ). It 
is important to note that the constant L defines directly 
the transmittance T of the superconducting slab. Indeed, 
according to Eq. ([T9|) . we have 



hi 



hi 



(22) 



The nonlinearity of Eq. (J16p leads to the multivalued 
dependence of the transmittance on the amplitude of the 
incident wave. In the following subsections, we analyze 
this unusual dependence, for both cases of negative {VL < 
1) and positive (Vl > 1) frequency detunings. 



B. Transmittance of a superconducting slab for 

CO < (Jjj 

We start with the case when the frequency of the in- 
cident wave is smaller than the Josephson plasma fre- 
quency. In this frequency range, the linear Josephson 
plasma waves cannot propagate in layered superconduc- 
tors. This corresponds to an exponentially small trans- 
mittance in the slab, due to the skin effect. However, the 



nonlinearity promotes the wave propagation because of 
the effective decrease of the Josephson plasma frequency. 

Solving Eq. (|16|) with the boundary condi- 
tions ([T5| . ([T^ . and (PH]) one can find the constant L and 
then calculate the transmittance using Eq. (|22p . Figure [2] 
demonstrates the numerically-calculated dependence of 
the transmittance T on the normalized amplitude ho 
of the incident wave. To analyze this dependence, we 
consider the spatial distribution of the gauge-invariant 
phase difference ip of the order parameter and the 
phase trajectories a' {a). We show these curves a' (a) in 
Fig. [31 An increase of the spatial coordinate C, [which 
is essentially z, as defined in Eq. (fT3)) ] from zero to 
S corresponds to moving along the phase trajectory 
a' {a). The point C = (i-S-, z = 0) corresponds to the 
starting point on the phase trajectory a' (a). According 
to Eq. ((20)) . all phase trajectories start from the points 
where a' = 0. Different trajectories in a' versus a can 
be characterized by the values of a(0) in these starting 
points. Each trajectory corresponds to some value of 
the normalized amplitude ho of the incident wave, and, 
according to Eqs. ([TS]), (HH), and 1^, the value of a(0) 
defines the constant L and the transmittance of the 
slab. 

The low-amplitude (quasi-linear) branch of the T{ho) 
dependence ranges within the interval < ho < (8/27)"'^/^ 
of the amplitudes of the incident waves. This branch is 
shown in Fig. [2] by the red solid curve close to the ab- 
scissa. For small amplitudes /iq ^ 1, we deal with a 
linear problem, when the phase difference ip and the elec- 
tromagnetic field in the superconductor can be found in 
the form of linear combinations of exponential functions 
of z. In this case, the transmittance T can be found 
asymptotically for small F, 



T{ho « 1) 



4r4 



sinh^[(5(l-K-2)]+4F4' 



F < 1. (23) 



This transmittance is very close to zero regardless of the 
frequency detuning (fi — 1). As we will see below, the 
"sinh" above, for w < loj, will become "sin" for u > ujj. 

The phase trajectories that correspond to the low- 
amplitude solutions occupy the region a < (8/3)^/^. 
For small ho, these trajectories are close to the point 
(a = 0, a' = 0) (as an example of such trajectory, see 
curve 6 in Fig. [3]). An increase of the amplitude ho leads 
to the growth of the length of the phase trajectory and 
this length tends to infinity when ho — > (8/27)^/^ (see 
curve 1 in Fig.[3|). 

The high-amplitude branches of the T{ho) dependence 
correspond to the solutions with a{C,) > 8^^^. Such 
branches arc shown in Fig. [2] by dotted, dashed, and blue 
solid curves for different values of the frequency detuning. 
The high-amplitude solutions describe nonlinear Joseph- 
son plasma waves that can propagate in the layered su- 
perconductor even for negative frequency detuning (for 
ri < 1). The corresponding phase trajectories are closed 
curves (see the closed curves in Fig.[3]for a > 8^^^). Note 
that the value of h is negative for a > 8^^^. For this case, 
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FIG. 2: (Color online) Dependence of the transmittance T 
[see Eq. (|22|) ] on the normalized amplitude ho of the incident 
magnetic field [see Eq. (|21l) ] for different values of the negative 
frequency detuning: (fi — 1) = — 5-10~^ (solid curve), — 5-10~ 
(dashed curve), —5 • 10"'^ (dotted curve). Arrows show the 
change of the transmittance when changing /iq. The numbers 
near the points on the solid T{ho) curve correspond to the 
same numbers of the phase trajectories a' (a) shown in Fig. [31 
The inset shows the enlarged region near the point 1. The 
values of the parameters are: S = 2, Xc — 4 ■ 10^'^ cm, Xab ~ 
2000 A, ujj/2Tr = 0.3 THz, and = 45°. 



we can consider h to be positive, but the phase of the 
incident wave must be shifted by tt. 

The oscillating character of the high-amplitude solu- 
tions results in much higher values of the transmittance, 
compared to the case of exponentially-small quasi-linear 
solutions. As seen in Fig. [51 the transmittance varies 
over a wide range, from nearly zero to one, depending 
on the amplitude ho of the incident magnetic field. It is 
important to note that the wavelengths of the nonlinear 
waves in the superconductor depend strongly on the inci- 
dent wave amplitude /iq. So, changing ho one can control 
the relation between the wavelength and the thickness of 
the slab. The transmittance is very sensitive to this re- 
lation, and one can attain total transparency of the slab 
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FIG. 3: (Color online) Phase trajectories a' {a) for the negative 
frequency detuning SI — 1 = —5 • 10~^. The numbers near 
the curves correspond to the same numbers of points in the 
T{ho) plot shown by the solid curve in Fig. [2] The movement 
along the phase trajectories in the directions shown by the 
arrows corresponds to the growth of the spatial coordinate 
1^ (proportional to z), from zero to 5, inside the slab. The 
solid lines show the portions of the phase trajectories that 
correspond to < (^ < S. The lower and upper surfaces of the 
slab correspond to the solid circles on the trajectories. The 
inset in the bottom left shows the enlarged region near the 
point (a = 0, a' = 0). The other parameters are the same as 
for Fig. [3 



choosing the optimal value h™'^^ of the amplitude ho . 

For high enough amplitudes ho, the sample thickness 
D is larger than the half-wavelength. In this case, the 
change of the coordinate C in the interval < ( < 6 cor- 
responds to the movement along a section of the phase 
trajectory loop (see the trajectories 2, 3, 4 in Fig. [3]). 
When decreasing ho-, the wavelength increases, the move- 
ment along the phase trajectory approaches the complete 




FIG. 4: (Color online) Spatial distribution of the amplitude h 
of the magnetic field inside the superconducting plate. Solid 
and dashed curves are plotted for the points 5 and 1 shown in 
Fig. [2l respectively. Point 1 corresponds to the low-amplitude 
branch of the T{ho) dependence, when the transmission coef- 
ficient is close to zero, and the amplitude h of the magnetic 
field near the lower interface {z = 0) of the slab is exponen- 
tially small. The transmittance T for the point 5 is equal to 
one, and the amplitudes of the fields near the upper and lower 
interfaces of the slab coincide. Here, the frequency detuning 
is (0,-1) — —5 • 10~^, and the other parameters are the same 
as for Fig. [2] 



loop, and the transmittance of the slab increases. Finally, 
for a specific value of ho, the wavelength becomes equal 
to the sample thickness, the phase trajectory draws a 
complete loop, and the transmittance becomes equal to 
one (see the trajectory 5 in Fig. [3] and point 1 in Fig.[2|). 
The amplitude dependence of the transmittance can be 
found asymptotically for small values of the parameter F 
and for not very thick slabs, S <1, 



Tiho) 



4F^ 



S^ V 72/10 



amplitude reaches the critical value (8/27)^/^ (point 1), 
further movement along this branch is impossible, and 
a jump to point 2, on the high-amplitude branch, oc- 
curs. A further increase in the amplitude ho results in a 
monotonic decrease of the transmittance. 

Afterwards, if the amplitude ho starts to decrease, then 
the T{ho) dependence does not follow the same track. 
Indeed, when the point 2 is passed, the transmittance 
continues to follow the high- amplitude branch. Decreas- 
ing the amplitude ho results in a further increase of the 
transmittance. When it becomes equal to one (point 5), 
it is not possible to continue the movement along the 
high-amplitude branch, and a return jump to the low- 
amplitude branch occurs. 

It should be noted that the jump from the low- 
amplitude branch (which corresponds to the exponen- 
tially small transmittance) to the high-amplitude one 
(with much higher transmittance) can be observed when 
changing the wave frequency w for the constant ampli- 
tude Ho of the incident wave. This jump occurs when the 
frequency detuning (1 — Vi) becomes equal to the thresh- 
old value 



(1 -f^ci) = - {^ckx:r^ 



2/3 



C. Mechanical analogy 

The problem discussed in this paper has a deep and 
very interesting mechanical analogy. Indeed, Eqs. (|16p 
and Eq. (J14p describe a motion of a particle with unite 
mass in a centrally symmetric potential. The amplitude 
/i(C) of the magnetic field, the phase 77(C), a-nd the coor- 
dinate C across the layers of the superconductor play the 
roles of the radial coordinate of the particle, its polar an- 
gle, and time, respectively. Moreover, the constant L in 
Eqs. (fT6|) and Eq. ([T4l) can be regarded as the conserved 
angular momentum of the particle. 

Integrating Eq. ((T6)) for the radial motion of the parti- 
cle, we obtain the energy conservation law for the parti- 
cle. 



(- 



(24) 



[h' 



U,ffih) = £, 



(25) 



In the case of total transparency {T ~ 1) of the slab, 
for h^^^ ~ 2^/^r(5, both the electric and magnetic fields 
take on the same values on the upper and lower surfaces 
of the slab. Thus, the amplitudes of the incident and 
transmitted waves are equal. The corresponding spatial 
distribution of the magnetic field is shown by the solid 
curve in Fig. |4l 

The nontrivial feature of the T[ho) dependence can be 
seen from its hysteretic behavior with jumps. Let the am- 
plitude ho of the incident wave increase from zero. In this 
case, the transmittance is close to zero, and the Tiho) 
dependence follows the low-amplitude branch shown by 
the red solid line near the abscissa in Fig. [51 When the 



with the effective potential 



a{h) dh. 



(26) 



The first term in Eq. ((25)) describes the kinetic energy of 
the radial motion of the particle, £ is the total energy of 
the particle. The first term in Eq. (|26p is the centrifugal 
energy and the last two terms represent the potential of 
the central field. 

The plot of the dependence Ucsih) is shown in Fig. [5] 
for the case of negative detuning (51 < 1). This de- 
pendence is three-valued and corresponds to the three 
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FIG. 5: (Color online) Dependence of the effective potential 
C/cff defined by Eq. H26|) on the radial coordinate h. The move- 
ment of the particle in this potential represents the mechanical 
analog for the spatial distribution of the amplitude h of the 
magnetic field in the superconductor. The main panel shows 
curves I, II, and III that correspond to the three branches of 
the a{h) dependence for the case of negative frequency detun- 
ing, n < 1. The inset shows the Ucnih) curve for the opposite 
case, Q, > 1, when only one branch of the a{h) dependence 
exists. The value of constant L is 0.1. 



branches of the function a{h) (see Eq. p3|) '). Thus, the 
multivaluedness of the dependence a{h) results in the 
multivaluedness of the cfFectivc potential Ucs{h) and, 
therefore, there exist several possibilities for the parti- 
cle motion. In terms of our electrodynamical problem, 
this means that several field distributions in the super- 
conductor can be realized for the same amplitude /iq of 
the incident wave. 

Curve I in Fig. [5] shows the potential that corresponds 
to the low-amplitude solutions of our electrodynamical 
problem. The motion of the particle (from right to left) 
in this potential is monotonic that corresponds to the 
monotonic decrease of the field deep into the supercon- 
ductor. According to Eq. (|20|) . the stop-point of the par- 
ticle {h' = 0) corresponds to the lower boundary of the 
superconductor. 

Since curve I is terminated in the point h = (32/27)^/^, 
it cannot define the particle motion for high enough h. 
In this case, the particle moves in the potential described 
by curve II in Fig. [S] This motion is finite and periodi- 
cal. It corresponds to the high-amplitude solutions of our 



electrodynamical problem. Curve III in Fig. [5] represents 
a branch of the Ucs{h) dependence which cannot be re- 
alized when changing the amplitude /iq of the incident 
wave. 



D. Transmittance of a superconducting slab for 

Now we study the transmittance of a slab of layered 
superconductors for waves with frequencies higher than 
the Josephson-plasma frequency, 11 > 1. Contrary to 
the case 17 < 1, even linear Josephson plasma waves can 
propagate in the layered superconductor. Therefore, the 
transmittance is not exponentially small and can vary 
over a wide range, depending on the relation between 
the wavelength and the thickness of the slab: 



T{hQ « 1) 



4r4 



sin-'[(S(l-K-2)]+4r4' 



r<i. (27) 



Note that the "sinh" in Eq. (P^ . for cj < wj, has now 
been replaced by a "sin" in Eq. (p7)) . for oj > ojj. 

In the nonlinear case, changing the amplitude ho, one 
can control the relation between the wavelength and the 
thickness of the slab and, thus, the transmittance is tun- 
able by the amplitude of the incident wave. Figure [6] 
shows the T{ho) dependences for different positive fre- 
quency detunings. 

The analysis based on Eqs. ([T5|) (with choice of the 
sign "-"), dUl), IHl), dUl), ^, and g^ shows that 
the dependence T{ho) is reversible when the frequency 
detuning is larger than some threshold value, which is 
defined by the asymptotic equation 



(f^thr - 1) 



fDy/I sin 9 



r<i. 



(28) 



V V2TTXab 

An example of such a reversible T{ho) dependence is pre- 
sented by the dotted curve in Fig. [6] 

The hysteresis in the T(ho) dependence appears for 
frequency detunings smaller than the threshold value: 

n < flthr- (29) 

In this case, the transmittance can reach the value one 
when the incident wave amplitude ho is first increased 
and then decreased. Namely, the incident wave ampli- 
tude ho is decreased after it increases and a jump of 
T(/io) occurs from the low-amplitude branch to the high- 
amplitude one (see the solid curve and the inset in Fig.[n|). 
One can derive the asymptotic equation for the optimal 
value h^'^^ of ho when the superconducting slab becomes 
totally transparent, 

^^r6^ (30) 



u max 
''0 



4/2 



where 



/ = 



dx 



VT 



,4/3 



~4^l2'4 



1.7972 
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FIG. 6: (Color online) Dependence of the transmittance T on 
the normalized amplitude ho of the incident wave for different 
positive values of the frequency detuning: — 1 = 5- 10"'^, 
or(5/7r = 1.2 (dotted curve); f7- 1 = 4.5 ■ 10-^ or S/-k = 1.25 
(solid curve); 0,-1 = 1.65-10"^, or S/-k = 2.1 (inset). Arrows 
show the change of the transmittance when changing ho . The 
sample thickness is _D = 4.3 ■ 10""" cm, and other parameters 
are the same as in Fig. (2] 
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FIG. 7: (Color online) Solution of the inverse problem: de- 
pendences of the amplitude ho of the incident wave and trans- 
mittance T = h^/ho on the amplitude hr of the transmitted 
wave. The values of the parameters and the numbers near 
the indicated points are the same as in the main panel in 
Fig. [S] The dependences plotted by dotted curves are mono- 
tonic, leading to the single- valued dependence T(ho) (the dot- 
ted curve in Fig. [S]). The solid-and-dashed curves show the 
nonmonotonic behavior that results in the multivalued depen- 
dence T{ho) (the solid curve in the main panel in Fig. [6|. 



and B(a;, y) is the Eulcr integral of the first kind. 



E. Mechanical analogy revisited 

Returning to the mechanical analogy described in the 
previous subsection, we note that, in the case of positive- 
frequency detuning, the dependence of the potential C/eff 
on the radial coordinate h of the particle is single- valued 
(see the inset in Fig. [5]). This is because the dependence 
a{h) in Eq. P^ is single-valued in this case. Neverthe- 
less, the dependence r(/io) can be multivalued even for 
il > 1 (see Fig. [5]). This feature seems to be paradoxical. 
Indeed, the particle motion is completely defined for any 
initial conditions. However, an assignment of the value 
of ha in relations ([TSl). (|19p . and ([20)) does not mean an 
imposition of definite initial conditions for the particle 
motion. To illustrate this nontrivial feature of the elec- 
tromagnetic wave transmission through a slab of layered 



superconductor, let us now consider the inverse problem. 
We wish to find the amplitude /iq of the incident wave 
that is necessary to obtain a given value Ht of the trans- 
mitted wave. According to Eq. ([T^ . the value of hx de- 
fines unambiguously the angular momentum L ~ h'^/T 
of the particle. On the basis of the motion equation (fTB]) 
and the boundary condition Eq. p8)) . we find that the 
dependence hQ^hx) should be single- valued. Correspond- 
ingly, the dependence of the transmittance 



T 



h^T/K 



on the amplitude Ht of the transmitted wave is also 
single- valued (see Fig. [7]). However, this dependence is 
nonmonotonic if the condition Eq. (|29p is satisfied. As 
a result, the dependence Tih^) appears to be multiple- 
valued. 



IV. CONCLUSION 

In this paper we predict a new nonlinear phenomenon 
in layered superconductors. We show that the reflectiv- 
ity and transmittance of a superconducting slab are very 
sensitive to the amplitude of the incident wave due to the 
nonlinearity of the Josephson relation for the c-axis cur- 
rent. As a result, the reflectivity and transmittance vary 
over a wide range, from nearly zero to one (if neglect- 
ing the small dissipation), when changing the amplitude 
of the incident electromagnetic wave. A remarkable fea- 
ture of this phenomenon is the hysteretic behavior of the 
r(ft,o) dependence. It is important to note that, for fre- 
quencies close to the Josephson plasma frequency, the 
tunable transmittance can vary from zero to one even in 



the case of weak nonlinearity, when the interlayer phase 
difference ip is small, i^ ^ 1. 
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